Dynamical entanglement and chaos: the case of Rydberg molecules 



O 
O 

(N 

G 

1—5 

00 



M. Lombardi and A. Matzkin 
Laboratoire de Spectrometrie physique (CNRS Unite 5588), 
Universite Joseph-Fourier Grenoble- 1, BP 87, 38402 Saint-Martin d'Heres, France 

A Rydberg molecule is composed of an outer electron that collides on the residual ionic core. 
Typical states of Rydberg molecules display entanglement between the outer electron and the core. 
In this work we quantify the average entanglement of molecular eigenstates and further investigate 
the time evolution of entanglement production from initially unentangled states. The results are 
contrasted with the underlying classical dynamics, obtained from the semiclassical limit of the core- 
electron collision. Our findings indicate that entanglement is not simply correlated with the degree 
of classical chaos, but rather depends on the specific phase-space features that give rise to inelastic 
scattering. Hence mixed phase-space or even regular classical dynamics can be associated with high 
entanglement generation. 
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I. INTRODUCTION 

Simply stated, a Rydberg molecule is composed of a 
highly excited electron orbiting around a compact ionic 
molecular core, containing the nuclei and the tightly 
bound other electrons. Most of the time, the outer 
electron is very far from the core. It is spatially well- 
separated from it, and the core and electron dynamics 
are uncoupled, the core's dynamics being characterized 
by its rotational energy. However, the outer electron pe- 
riodically scatters on the molecular core. In quantum- 
mechanical terms the electron and core dynamics get cou- 
pled and the collision induces phase-shifts in the wave- 
function of the outer electron, known as quantum defects 
0, Q . This process has a well-defined classical counter- 
part Q: the electron is kicked by the core, resulting in 
a change of the electron's angular momentum relative to 
the core by an angle which depends on the quantum de- 
fects. During the kick the core and the outer electron may 
exchange energy, so that in general the rotational state 
of the core has changed after the collision. Although a 
molecule is intrinsically a quantum object, several prop- 
erties, such as the stroboscopic effect seen in laser excita- 
tion or the statistics of the energy levels have been shown 
to depend on the underlying classical dynamics 0, Q . 

There is however a distinctive quantum feature that 
has no classical counterpart, which is readily touched 
upon by observing that a typical molecular state is de- 
scribed by quantum mechanics as a superposition of the 
different core states available to the molecule, each state 
being associated with the corresponding outer electron. 
This superposition is the result of the electron-core en- 
tanglement produced by the collision. 

The main goal of this paper is to investigate the en- 
tanglement dynamics and in particular its dependence on 
the underlying classical regime. Indeed there has been 
a growing interest in recent years to correlate the en- 
tanglement production of a quantum system with the 
dynamics of the corresponding classical system. For ex- 
ample the time-dependent entanglement production was 
investigated on the N atom Jaynes-Cummings model (a 



single mode field interacting with a 2-level spin), where 
the initial product-state wave-packet was chosen to lie 
in different regions of classical phase-space ttl. Similar- 
studies were undertaken for kicked tops 0, 0, ■ The 
initial claims 0.0 by which entanglement would system- 
atically increase with chaos were revised in subsequent 
works. One of the problems was to unequivocally define 
the classical counterpart for such systems, which is far 
from being straightforward for the spin-boson case and 
has been discussed for the coupled kicks system. A semi- 
classical approach based on the average properties of 
phase-space concluded that both the coupling strength 
(in the form of classical correlators) and the global clas- 
sical dynamical regime were important in understanding 
the entanglement rate. However the relevance of univer- 
sal relations relating the generation of entanglement to 
the dynamics of the corresponding classical system are 
still being debated: the form of the initial quantum state 
is known to play a role 0, ITTl IT^ . and from a more 
general viewpoint realistic systems usually display more 
complex dynamics than simple systems that present a 
uniform behavior over all points of phase-space. In this 
respect the choice of Rydberg molecules to investigate 
entanglement generation may be fruitful: these real sys- 
tems are theoretically well described (by quantum defect 
theory), the semiclassical limit is readily obtained, and 
the classical dynamics is sufficiently simple to be well un- 
derstood (at least qualitatively) without being too sim- 
plistic. 

This work is organized as follows. In Sec. II we briefly 
recall some basics concerning the quantum theory of sim- 
ple Rydberg molecules and describe the classical counter- 
part of such systems, insisting on the relevant Poincare 
surfaces of section. In Sec. Ill we determine the degree of 
the electron-core entanglement, inferred from the linear 
entropy of the reduced density matrix for the outer elec- 
tron. We will calculate some simple statistics on groups 
of eigenstates corresponding to different classical regimes. 
We will then follow the time-evolution of the entangle- 
ment from an initially unentangled product state. This 
will be done for different initial states and various dynam- 
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ical regimes. The results will be discussed in Sec. IV. We 
will sec that the classical dynamics is reflected in the en- 
tanglement generation; but the global dynamical regime 
is less important than the specific classical effects that 
are quantum mechanically translated into superpositions. 
Our closing comments will be given in Sec. V. 

II. RYDBERG MOLECULES - QUANTUM AND 
CLASSICAL 

A. Quantum phase-shifts 

A Rydberg molecule is composed of two elements: a 
highly excited electron - the Rydberg electron, on the one 
hand, and the ionic molecular core, positively charged, 
containing the nuclei and the tightly bound other elec- 
trons on the other hand. The Rydberg electron is usually 
far from the core and only senses the long-range Coulomb 
field produced by the core, irrespective of the complex in- 
teractions involving the core particles. The total energy 
of the molecule E is consequently partitioned between 
the energies of the core and of the Rydberg electron. 

For definiteness, we will take a diatomic molecule for 
which vibration can be neglected: rotation is then the 
only motion available to the core. Since the molecule is 
isolated, the total angular momentum J and its projec- 
tion M on an axis fixed in the laboratory are conserved. 
We have 

J = N + L, (1) 

i.e. J results from the addition of the angular momenta 
of the core, N, and of the Rydberg electron L. We shall 
assume that L is conserved, as is often the case, and the 
standard addition of angular momenta gives 

\J-L\ <N< \J + L\. (2) 

Each state of the core is therefore labeled by \N). Recall 
that it also follows from the addition of angular momenta 
that whereas N and L are well-defined, their projections 
Mjv and m on the reference axis are not, since only the 
total projection M = Mjv + m is well-defined. Therefore 
the notation \N) also contains the angular state of the 
Rydberg electron, via the geometrical angular momenta 
couplings. As is well known, the energy of the core E + 
depends on N through a rotational constant B r , 

E+=B r N(N + l), (3) 

and the total energy of the molecule is thus 

E = B r N(N +l) + e N (4) 

where e./v is the energy and vn the effective quantum 
number of the Rydberg electron, 

-1 

€n = ~o77T (5) 

N 



(atomic units are used throughout). Note that for a given 
value of E, e implicitly depends on the state of the core: 
the electron is more or less excited depending on whether 
the core has a large or small rotation number N. Here we 
will only deal with bound states (i.e., E is below the low- 
est ionization threshold). The wavefunction correspond- 
ing to the Rydberg electron in the Coulomb field of a 
core in the rotational state \N) is thus 

ct> N (E,r)=\N)f L (E-E+,r), (6) 

r being the radial coordinate and /l the Coulomb func- 
tion regular at the origin. 

Now when the electron significantly approaches the 
core, the interaction between the core particles and the 
electron cannot be neglected. These interactions are em- 
bodied in the short-range potential V, that is negligible 
beyond the core radius ro- Thus the total Hamiltonian 
of the molecule is 

H = H + V, (7) 

with the functions <f>N being eigenstates of Hq. From 
scattering theory, it is well-known that the wavefunctions 
for H are obtained from 

4> N {E, r) = <j> N {E, r)-J29L(E- E+„r) \N') K N , N , 

N' 

(8) 

where <7l is the Coulomb function irregular at the origin 
and K is the scattering matrix; Kn>n gives the transition 
probability amplitude between states 4>n and 4>N' during 
the collision. The eigenfunctions of the total Hamiltonian 
H are obtained by the superposition 

iP(E,r) =Y,Z N (E)^ N (E,r). (9) 

N 

The coefficients Z^(E) are obtained from quantum de- 
fect theory Q by imposing the appropriate boundary 
conditions at infinity, yielding also the discrete eigenval- 
ues E. We thus see from Eqs. J8J and JJjJ that a generic 
wavefunction for a molecular Rydberg state involves a su- 
perposition of core states with different rotational num- 
bers. This is caused by the short-range potential V that 
transforms the product state </>jv [Eq. JHJ] into the en- 
tangled state tp- I n practical computations, Eq. JjJ) i s 
of awkward use because taken individually each of the 
functions ipN(E,r) diverges radially. Eq. is therefore 
rewritten as 

iP(E,r)=J2 B N(E)F L (e N ,r)\N) (10) 

N 

where the FL(eN,r) are the effective channel functions, 

F L (e N ,r) = sm/3(e N )f L (e N ,r) - cos P(e N )g L (e N , r), 

(11) 

which by construction converge as r — > oo. /3(ejv) = 
7r(^jv — L) is precisely the total phase accumulated at 



FIG. 1: Molecular reference frame. The core axis is along OZ 
and its angular momentum N is along OX 



r — > oo. The coefficients Bjq are obtained from the 
by matching Eqs. © and (jTU|) . 

An additional subtlety arises from the use of frame 
transformations: the wavefunctions given above were de- 
scribed in the laboratory frame. However when the elec- 
tron collides on the core, a description in the molecu- 
lar frame, attached to the core rotation, is more appro- 
priate, because in the core region the Rydberg electron 
senses the cylindrical field aligned along the molecular 
axis. Thus only the projection of L on that axis, tradi- 
tionally denoted by A, matters when describing the col- 
lision: the phase-shifts induced by the collision on the 
Rydberg electron's wavefunction are known as quantum 
defects and denoted by fj,\. The collision matrix K ap- 
pearing in Eq. © is obtained by expressing the phase- 
shifts ha in the laboratory frame 0. 



B. Classical kicks 

The classical counterpart of the quantum model in- 
troduced above is the following When the Rydberg 
electron is far from the core, it follows a pure Coulomb 
(Kepler) orbit with an angular momentum L fixed in 
space. Meanwhile the molecular core rotates freely with 
a rotational energy N(N +1)/2I depending on the core's 
angular momentum N and moment of inertia / = l/2B r . 
Seen in the molecular frame, L precesses around N, i.e. 
it turns around the N axis with a constant angle (see 
Fig. Now when the outer electron approaches the 
core, it gets kicked by the molecular axis. This kick re- 
sults in a change in the direction of L. N adjusts accord- 
ingly, since the total angular momentum J is conserved. 
To visualize the effects of the kick, it is therefore sufficient 
to follow the evolution of L in the molecular frame. 

Since during the collision the Rydberg electron feels 
the cylindrical field due to the molecular axis, 8 cannot 
change and thus only the angle ip varies (see Fig. [Q. 
This variation, denoted 6<p, is the deflection angle of the 
plane of the classical Kepler orbit. The relation between 
the classical deflection angle and the quantum scattering 



FIG. 2: Two limiting Poincare surfaces of sections, (a) no 
coupling k — (b) large coupling k — 10 (in the generic case, 
see Fig. 0] below) . The shaded areas correspond to quantum 
states with minimal ( J — L = 40), and mean (J = 50) values 
of N (see Sec. III). Their widths correspond to AN = 1. 
The mean N — J case is not on the equatorial OYZ plane, 
as can be understood by squaring J = L + N (where N is 
"horizontal", parallel to OX) when the lengths of N and J 
are equal. 



phase-shifts is well known from the semiclassical approx- 
imation to the scattering amplitude Here it takes 
the form 



Sip = 2tt 



OA ' 



(12) 



i.e. the strength of the kick, reflected in the amplitude 
of the deflection angle, is the classical counterpart of the 
dependence of the phase-shifts on A. The precise form 
of this dependence depends on the particular molecule at 
hand. However, for a typical molecule, the dependence 
can be taken in the form 



MA = Mo - k 



A 2 



4ttL 



(13) 



Classically the coupling parameter k gives the strength 
of the kick since it follows from Eqs. I|12f> and (|13fl that 



\5(p\ = k— = k cos ( 
L 



(14) 



To visualize the classical dynamics, a Poincare surface 
of section is obtained by plotting the position of L after 
each kick. This is most naturally done in the molecu- 
lar frame where the position of L is given by the polar 
angles (0,ip), as shown in Fig. ^ with ip being canon- 
ically conjugate to L z . This is the molecular reference 
frame used in quantum mechanics. However, as in our 
previous works |3, 0, [w| , we prefer to add an extra ro- 
tation around OZ to bring the OX axis along N. This 
last rotation is not canonical but the classical motion is 
seen more naturally in this frame (see the discussion in 
Ref. H3)- The position of L is thus plotted after each 
kick, when the electron comes out of the core. Of course, 
since J = L + N is conserved, following the position of L 
is tantamount to knowing the fate of N. Two extreme 
examples of Poincare sections are given in Fig. [2J In 
the first case, Fig. E[a), the positions of L follow circles 
around N (the X axis) giving an overall regular surface 
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of section. This means that on a given circle L is fixed in 
space, as is always J, and thus that N = |J — L| is con- 
stant. In the second case, Fig.^b), the surface of section 
is clearly characteristic of a chaotic phase space and no 
structure arises by following the successive positions of L 
after each kick. 



III. ENTANGLEMENT DYNAMICS 

A. General Remarks 

Generically, in a Rydberg molecule the core and the 
outer electron are entangled; in the situation examined 
in this work, the core is in a superposition of different 
rotational states. Each rotational state is defined by a 
given value of N, and to each core state is associated 
an outer electron with an energy e(N) given by Eq. 
Hence the outer electron is in a superposition of different 
channels which arc distinguished in the present model by 
a different energy, depending on AT. Classically of course 
there is no superposition in N: the rotational energy of 
the core, and thus the energy of the outer electron, is 
at each instant unique and well defined. A change in 
N can only be the dynamical result of an inelastic kick. 
Note that classically as well as quantum mechanically, 
the angular momenta L and N are coupled, since J is 
conserved. As is well-known, the composition of angu- 
lar momenta in quantum mechanics results in couplings 
due to the fact the projections of the angular momenta 
vectors cannot be simultaneously defined. This entangle- 
ment of geometrical nature (as in EPR pairs) should not 
be confused with the dynamical entanglement generated 
by the potential interaction between the outer electron 
and the core. Geometrical entanglements due to angular 
momenta coupling do not play any role in this work, as 
we are only interested in the dynamical one which only 
depends on |N|. This implies that a product state is given 
by the channel functions </>jv(-E, r), defined by Eq. (Jj), i.e. 
only the radial coordinate of the outer electron is separa- 
ble from the core whereas its orbital angular momentum 
is necessarily geometrically coupled to the core's angular 
momentum N (since J is conserved and Hq has spher- 
ical symmetry). Hence when referring to partial traces 
on the outer electron we will mean a trace over its sole 
radial coordinate and conversely a partial trace over the 
core includes tracing over the angular coordinates of the 
electron. 

To quantify entanglement we will determine the linear 
entropy 5*2 associated with the reduced density matrix 
p e descibing the outer electron, 

Pe = Tr cP = ^2(N\p\N), (15) 

N 



where p = \tp) (tp\ is the density matrix of the system and 
Tr c (resp. Tr e ) refers to averaging over the core (resp. 
outer electron) degrees of freedom. The reduced linear 
entropy is then defined by 

S 2 = l-Tr e p 2 e . (16) 

Strictly speaking, 5*2 measures the degree of mixedness: 
S2 vanishes for a pure state and is maximum for a uni- 
formly mixed state. However when S2 is associated with 
improper mixtures in bipartite systems, it reflects the de- 
gree of entanglement |15| and the linear entropy or equiv- 
alent quantities such as the pur ity have routinely been 
employed as such 0, 0, 9, 10, 11, 161. We will undertake 
two different studies. First we will investigate the degree 
of entanglement on stationary states and its dependence 
on the collision phase-shifts, whose classical counterpart 
gives rise to different dynamical regimes. This involves 
the computation of simple statistics of S2 (E) in an energy 
range for which the classical dynamics does not vary. We 
will then investigate the time evolution S2 (t) from an ini- 
tially (at t = 0) product state. This involves the determi- 
nation of wavepacket dynamics. The initial wavepacket 
can be made to lie and then evolve in zones corresponding 
to different classical dynamics. 

We will obtain numerical results for the following 
choice of parameters: J = 50, L — 10 yielding by Eq. 
21 values for N. Since Kronig's parity, i.e. parity by 
reflection at a plane through the internuclear axis OZ 
(which for fixed L and J depends on the parity of N), is 
conserved, entanglement only takes place between states 
of the same Kronig's parity. States of (+) and (— ) Kro- 
nig's parity, called positive and negative states in stan- 
dard spectroscopic notation 1 behave in the same way, 
and we will restrict our study to (+) parity states; then N 
can only take even values, so that a typical state contains 
superpositions involving up to 11 values of N. These an- 
gular momenta numbers are considerably higher than for 
typical diatomic molecules (they would better correspond 
to models of large molecular compound Rydberg states). 
However higher quantum numbers allow a finer compar- 
ison between quantum and classical dynamics, without 
qualitatively affecting the correspondence between them: 
see |l8| , where partial results with angular quantum num- 
bers typical of diatomic molecules such as Na2 were ob- 
tained 12311. 



B. Degree of entanglement of eigenstates 

The eigenstates tj){E) of a Rydberg molecule are given 
by [Eqs. (fTTTITll l 
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i>(E) = B N (E) (sin/3(£ - E+)f L (E - E+,r) - cos 0(E - E+)g L (E - E+,r)) \N) . 

I 



(17) 
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FIG. 3: Average linear entropy for generic (red boxes, color 
online) and resonant T e = To states (black boxes), shown for 
different values of k. The rms is also shown as red dashed 
(generic) and solid (resonant) error bars. 






FIG. 4: Poincare surfaces of section for different values of the 
coupling constant k. Top: generic; bottom: resonant cases. 
(a,d) k = 0.25. (b,e) k = 1 (c,f) k = 10. 



To quantify the degree of entanglement of a given eigen- 
state, we compute 82(E). Variations from individual 
eigenstates are smoothed out by calculating simple statis- 
tics for the bunch of eigenstates sitting in an interval AE. 
The requirement on AE is that the classical dynamics 
does not change appreciably within this interval. Several 
computations are performed for different values of the 
coupling constant k. Each value of k corresponds to a 
different collision matrix K via Eq. (|13|l . 

The results are shown in Fig. [3] The average and rms 
of 52 are given for different values of k. Fig. shows 
both the results for a 'generic' situation which would be 
obtained for an arbitrary choice of E and the results for 
resonant eigenstates: in the latter case E is chosen such 
that the period of the Rydberg electron T e is an integer 
multiple of half the period of the core T c . Classically, this 
corresponds to a situation in which the electron sees the 
core in the same position on its return as when it left the 
core region. Resonances affect the classical dynamics, 
essentially by retarding the appearance of chaos. This 
is portrayed in Fig. 0] which shows Poincare surface of 
sections for several values of the coupling k in both the 
generic and the resonant cases. In the former case, chaos 
appears even for a small value of fc, whereas in the latter 
configuration, chaos becomes significant for larger values 
the coupling, and even for such a large value as k = 10 
an island of regularity around the fixed point on the Z 
axis is still visible. 

Figure|21indicatcs that on average the degree of entan- 
glement, as measured by < £2 > increases with k. How- 
ever this does not mean that entanglement is correlated 
with classical chaos. Thus for example for k = 0.25 the 
average entanglement is significantly higher in the reso- 
nant case than in the generic case, despite phase space 
being slightly more regular (see Figs. 01a) and (d)). For 



k = 1 < S2 > has the same value in the generic and res- 
onant cases although the dynamical regimes, as reflected 
in the surfaces of section Fig.^Jb) and (e), are quite dif- 
ferent. 

We further illustrate the effect of a resonance both on 
the classical dynamics and on the linear entropy for the 
case k = 0.5. Indeed, as the energy is appreciably moved 
away from the exact resonance energy, the classical dy- 
namics accordingly changes, going back towards a generic 
situation. This is shown in Fig. (Top) : the Poincare 
surface of section at the center is plotted at the exact 
resonance energy. However as the energy changes appre- 
ciably, the periods of the Kepler orbit and the core rota- 
tion are significantly altered, suppressing the resonance. 
Classically, the structure of phase-space is modified, as 
can be directly seen on the surfaces of section plotted 
at both ends of the energy range. At all energies phase- 
space is regular, but the separatrices characterizing the 
resonance give way to island chains and curves organized 
around OX. Fig. [5] (Bottom) shows the linear entropy 
for the individual states lying within this energy range. 
It may be seen that as we move away from resonance the 
behavior of the linear entropy drops dramatically: most 
states show a lower degree of entanglement. The degree 
of entanglement clearly appears to be correlated with the 
changes in phase-space induced by the resonance. The 
findings presented in Figs. |21 and [5] will be discussed in 
Sec. IV, but we may note that these results indicate that 
entanglement is sensitive to the details of classical phase 
space, not only to the global dynamical regime. Note 
that the statistics for k = 0.5 shown in Fig. [3] were done 
from the 200 individual states lying within the dashed 
lines in Fig [SJ 
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FIG. 5: Classical phase-space and linear entropy in the 
k — 0.5 resonant situation. (Top) Poincare surfaces of sec- 
tions at the left, middle, and right ends of the energy interval 
shown at the bottom. (Bottom) Linear entropy of the indi- 
vidual eigenstates. The black dots between the dashed lines 
correspond to the states that entered the statistics shown in 
Fig. [3 Each eigenstate of energy E is labeled by v(E), the 
principal quantum number in the N = J channel. 



To grasp the relationship between quantum entangle- 
ment and the classical dynamics for individual eigen- 
states, we project the wavefunctions in mock phase-space 
on the surface of section by drawing the Husimi plots for 
two eigenstates having different values of We give in 
Fig. El the Husimi plot along with the 11 Bm coefficients 
of two particular wavefunctions lying very near the center 
of the resonance shown in Fig. [5] The first one is (acci- 
dentally) a nearly pure N — 40 wavefunction. Its Husimi 
plot shows that it is quantized near the +OX axis, as 
expected from the gray zone displayed in Fig.|3 Being a 
nearly pure product, it has thus a very low linear entropy, 
S2 = 0.03658. The second wavefunction is (also acciden- 
tally) quantized nearly on the +OZ axis. Its phase space 
extension is nearly the same, but the decomposition on 
the Bn basis spans more values of N, as shown in the 
lower part of the figure. This is understood by consid- 
ering the overlap of its Husimi plot with OA-centered 
circles, each such circle corresponding to a value of N. 
Its linear entropy is thus correspondingly much higher: 
S 2 = 0.6881. 



C. Dynamical evolution of entanglement 

We consider the time dependence of the generation of 
entanglement from an initial product state. Assume the 
system has been prepared so that at t = the core has 
a well defined rotational state \Nq) whereas the outer 
electron is radially localized at the outer turning point 
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FIG. 6: Husimi plots (Top) and Bn(E) coefficients (Bottom) 
for two particular wavefunctions in the k = 0.5 resonant sit- 
uation near the center of the interval plotted in Fig. |S] Left 
level, with u(E) = 315.4233, is nearly quantized on the min- 
imum value of N — 40 as shown by the Bm distribution. 
Compare it's Husimi plot with shaded areas in Fig. It has 
S 2 (E) = 0.03658. The right level, with v E = 315.6118, quan- 
tized nearby the +OZ axis, spans a greater number of Bn 
values. It's Husimi plot overlaps with a greater number of TV 
circles indicated in Fig.H It has S 2 {E) = 0.6881. 



of the Kepler orbit, several thousand atomic units away 
from the core. The wavepacket attracted by the Coulomb 
interaction moves towards the core and collides on it at 
t w T e /2. The collision results in entanglement, since the 
outgoing waves of the Rydbcrg electron are in a state of 
superposition, each outgoing channel being attached to a 
core in a different quantum state. Subsequent collisions 
with the core result in further entanglement, whereas the 
spreading of the radial wavepacket quickly results in a 
continuous core-electron interaction. 
We take the initial state to be 

^{t = 0,r) = F loc {r » rtp) ® |JV ) (18) 

where the radial function is 

fioc(r) =J2^ l(n ~ no)/2An]2R nL(r) (19) 

n 

with an appropriate normalization factor. i? n i(r) are 
the standard radial functions of the hydrogen atom (reg- 
ular Coulomb functions) and the Gaussian form of the 
coefficients are known to ensure localization 19]. tiq is 
chosen so that the central component of the wavepacket 
matches the energy of the corresponding classical regime 
under study. At later times the wavefunction is given by 

M r ) = E E B N {E)e- iEt F L {E - E+,r) \N) , (20) 

E N 

where 

B N {E) = B N (E)B No {E) (F L (E - E+)\ F loc ) (21) 
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where the coefficients Bn(E) and effective channel radial 
functions F L (E-E+) were given above (cf Eqs. (|10() and 
ijnj). The radial overlap {Fl(€n)\ Fi oc ) is determined 
analytically as a particular instance of the scalar product 
(F L (e)\ F L (e')) given by [U 



(F L (e)\ F L (e')) = 



sin7r(i^ — v') 



tt(v — i> r ) 

multiplied by the relevant normalization factors 



(22) 



In principle, the computation of the linear entropy as- 
sociated with the reduced density matrix is straightfor- 
ward. From 



P (t) = m))(m\ 



(23) 



we obtain the purity Tr e pl(t) and the reduced linear en- 
tropy Sa(t) as 



TWeW = / {r\ P l\r)r 2 dr 



Tr e pl(t) = ]T 



NN' 



£ e-^ E - E '^B N {E)B N -{E') (F L {E' - E+, \ F L (E - E+)) 



EE' 



(24) 
(25) 




t/T e 

FIG. 7: Short time variation of the linear entropy generated 
from the initial product state Fi oc (r ~ r tp ) ® | A^o = 40) for 
different collision matrices labeled by the value of the coupling 
strength k: from top to bottom k — 10 (color online, red), 1 
(black) and 0.25 (blue), t is given in units of the Kepler period 
T e . (a) The parameters are chosen so that the corresponding 
classical dynamics falls in the generic phase-space case, (b) 
Same as (a) for the resonant case T e = To. In both (a) and 
(b) the k = 0.25 curve is multiplied by a factor 5. 

The radial closure relation needs to be introduced in 
Eq. Ij24(l given that the effective radial functions that 
play the role of the basis are overcomplete. 



1. Short-time evolution 

We examine first the short time-evolution for three dif- 
ferent couplings (kick strength) - k = 0.25, k — 1 and 
k = 10 - and two different initial states distinguished by 
the value of Nq: Nq = J — L (that is the minimal value 




2 4 6 8 10 12 

t/T e 



FIG. 8: Same as Fig. |7| but with the initial product state 
given by F; oc (r ss r tp ) ® \No — 50) (the k = 0.25 curve is not 
multiplied by a factor 5). 



Nq can take, cf Eq. J2Jl) and Nq = J. The correspond- 
ing zones in the surfaces of section are shaded in Fig. 
Recall that on the surface of section a fixed value of A*o 
corresponds to a line circle around the X axis, whereas 
a quantum state in mock phase-space projected on the 
surface of section has a certain width, as seen in Fig. 
The width of the initial state can be roughly estimated 
by plotting the zone going from N — 0.5 to N + 0.5, cut- 
ting the sphere approximately into the number of values 
N can take (here 11). 

Fig. Eta) and (b) show the linear entropy as a function 
of time (in units of the Kepler period) for the generic 
and resonant cases when the initial state l|18|) is chosen 
with No = J — L = 40. Figure [S] gives the linear entropy 
when the initial state is taken with iVo = 50. In Fig. 
we observe that in both the generic and resonant cases 
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S2(t) = 1 — Tr e pg(i) increases with k: entanglement is 
produced more rapidly and saturates at a higher value. 
Comparing with the classical dynamics, we see that in 
the generic case, the rise in entanglement generation ac- 
companies the classical transition to chaos. This remains 
true to a certain extent in the resonant mixed 
phase-space turns progressively chaotic. 

However when the initial rotational state lies near the 
center of the surface of section (Fig. 01, the linear entropy 
takes very high values irrespective of the classical dynam- 
ical regime. This is particularly spectacular for k = 0.25, 
which jumps from negligible values in Fig. to crossing 
the k = 10 curve of the linear entropy in Fig. |^b). For 
the first few periods this increase of the linear entropy 
takes place in steps, reflecting the collision of the radi- 
ally localized electron wavepacket with the core at each 
half-integer value of t/T e (the wavepacket spreads radi- 
ally after a few periods). These findings will be discussed 
in Sec. IV but we may note again that as found for the 
eigenstates there is no simple relation between the global 
classical dynamical regime and quantum entanglement 
generation. 



2. Long-time evolution 

In Figs. [7| and [SJ the linear entropy appears to satu- 
rate after a few periods. For longer times S2 (t) is plotted 
in Figs. |5] (when the initial state has No = 40) and 1101 
(No = 50). In most of the cases (all the k = 10 curves 
and the k = 1 curves except in Fig. GJb)) S^i) appears 
to vary randomly around some average value. However 
for small kicks (the k = 0.25 curves, but also the k = 1 
curve in Fig.^Jb)) the repetition of oscillatory structures 
is clearly visible to the eye. These repetitions are due to 
partial wavepacket revivals within each channel, which 
take place when the terms that control the spreading of 
the packet regain an approximate phase coherence. The 
revival times within each channel are determined in the 
semiclassical approximation by expanding the energies in 
the exponentials in Eq. 112511 as a function of the classical 
action. The core's revival time T£ ev cx B~ 1 is indepen- 
dent of the energy whereas TJ e " cx e~ 2 . A revival in the 
linear entropy will be visible provided 1 or 2 channels 
dominate in the overall contribution. 

For the k — 0.25 case in Fig. EJb), most of the prob- 
ability density stays in the original No = 40 channel, a 
small flux flowing to the N = 42 channel. This flux is 
responsible for the main part of the entanglement gener- 
ation. We have plotted in Fig.^Jthe correlation function 



where 



C(t) = (tM* = r e /2)| M)) 



\^{t)) = {N = A2\^{t)) 



(26) 



(27) 
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FIG. 9: Long time variation of the linear entropy for the case 
considered in Fig. (No = 40, (a) generic and (b) resonant 
cases). The inset in (b) blows up the vertical scale for the 
k = 0.25 curve to help visualize the oscillations, also visible 
for k = 1. 
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gives the electronic wavepacket in the N = 42 channel. It 
can be seen that near t — 100 T„ the correlation function 
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FIG. 10: Long time variation of the linear entropy for the case 
considered in Fig. |H| (iVo = 50, (a) generic and (b) resonant 
cases). The inset in (b) details S^i) in the region t = 80 to 

110 T e . 



oscillates dramatically: high peaks appear while at the 
same time the lowest values are near zero. This behavior 
translates into the entanglement rate, which shows the 
same strong oscillations as seen in the inset in Fig. [5] 
The revival time for the electron motion in this case is 
computed in the semiclassical approximation as T£ ev ~ 
105 T e . This value fits well with the periodic repetitions 
of these strong oscillations observed for S2 it) . 



IV. DISCUSSION 

If we compare the time evolution of the linear entropy 
portrayed in Fig. to those of Fig. |H1 the importance 
of the initial state is clear. If we further contrast these 
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FIG. 11: The partial autocorrelation C(t) defined by Eq. 12611 
(parameters corresponding to the k — 0.25 resonant case, ar- 
bitrary units) is plotted in an interval centered around the 
revival time TI ev ^ 105 T e . 



results with the Poincare surfaces of section of Fig.0]and 
the localization of the initial states shown in Fig. |3 we 
see that the generation of entanglement does not essen- 
tially depend on the global classical dynamics, but rather 
on the specific dynamics that leads to inelastic scattering. 
Classically, inelastic scattering means that the value of N 
changes during the collision; then on the surface of sec- 
tion two consecutive points cannot lie on circles around 
the X axis (since these circles precisely correspond to a 
fixed value of N) . Inelastic collisions do take place when 
the dynamic is chaotic: for very strong kicks (k = 10) 
two consecutive points on the Poincare sections are arbi- 
trarily separated on the sphere. Quantum mechanically 
inelastic scattering is translated into superpositions of 
states having a different value of N. The k = 10 curves 
in Figs.[7|^a) and^a) indeed reflect large and fast entan- 
glement generation. S^i) achieves its maximal value (of 
10/11) just after a couple of collisions. 

However inelastic scattering can also be induced by 
regular dynamics. Consider k = 0.25 in the resonant 
case. When the initial state encircles the X axis at the 
front of the sphere as happens with No = 40 (see the 
shaded region in Fig. |2J), it will be hardly modified by 
the classical dynamics. The lines in the surface of sec- 
tion (Fig. Eld)) also encircle the X axis, thereby con- 
serving No: there is essentially only elastic scattering. 
Quantum mechanically we expect little or no entangle- 
ment, as is observed in Fig. [7| But for No = 50 the 
initial distribution spans across the lines of regularity in 
the surface of section, which are organized around the 
elliptic fixed point on the Z axis. These lines thus de- 
fine torii in phase-space that cut across several values of 
N, meaning inelastic scattering and quantum-mechanical 
superpositions. Note that in the generic k = 0.25 situ- 
ation the linear entropy is lower than in the resonant 
case; classically the structure of phase space in the zone 
covered by the initial state is modified due to the ap- 
pearance of resonant islands (Fig. Ufa)), leading to some 
lines that roughly encircle the X axis, partially favoring 
elastic scattering. 

We thus see that the entanglement production reflects 
the classical dynamics in that the relevant parts of space- 



phase leading to inelastic collisions (yielding quantum 
superpositions) are explored with a large amplitude. To a 
large extent this phenomenon is unrelated to whether the 
dynamical regime is chaotic or regular. Classical chaos 
is not necessary to induce high entanglement generation. 
In our system chaos is sufficient, because it always leads 
to inelastic scattering. Hence eigenstates quantized in an 
underlying strongly chaotic phase-space will present on 
average a large amount of mixtures, as seen in Fig. E3 
These remarks are in line with similar conclusions |l0l 
HTl ll2T ] which contrarily to earlier results 0, |(J , do not 
attribute to chaos a higher entangling power. 

Let us mention that more precise investigations of the 
correlations between classical dynamics and dynamically 
induced entanglement should take into account more in- 
formation than what can be inferred from the Poincare 
surfaces of sections. For example in the generic k = 1 
case (Fig.^Jb)) diffusion in the chaotic sea takes place at 
a considerably lower rate than in the A; = 10 situation. 
In particular, the particle may be trapped for several pe- 
riods in certain regions of phase-space, diffusing slowly 
in the relevant regions of the surface of section. This 
type of phenomena has an influence on the quantization 
process, and if important will influence the entanglement 
production, as in the present case. We therefore expect 
that these system specific features, along with the role of 
the initial state relative to the precise structure of phase- 
space, may severely constrain the applicability to realistic 
systems of general 'universal' formulas ruling the entan- 
glement generation in chaotic and regular systems that 
have been recently obtained 0, . Semiclassical uni- 
versal formulas, based on the global average properties 
of classical phase-space, are important as they set the 
trend that is followed by a quantum system with sim- 
ple dynamics (like a uniform transition to chaos). How- 
ever individual features of the system are known to be 
important in the semiclassical description of diffractive 
effects induced by a coupling potential and need to be 
taken into account e.g. to describe the spectral statis- 
tics [2l|. It is thus not surprising to see that dynamical 
entanglement induced by a standard potential coupling 
is correlated with the local structure of phase-space, and 
not only with its global properties. 

Concerning the long-time behavior, it may be noted 
that the revivals appear well beyond the Heisenberg time, 
time after which quantum phenomena having no classical 
counterpart become prominent. In some cases, for exam- 
ple when the number of channels is small, the oscillations 
due to wavepacket revivals induce very large variations of 
the entanglement rate totally unrelated to the behavior 
at short-times. Hence the linear entropy of a state which 
initially only showed a slow and small amount of entan- 
glement generation can rise above the linear entropy of a 
state that initially displayed high and fast entanglement 
production. 

We finally point out that entanglement in Rydberg 
molecules is routinely detected experimentally, given that 
the consequences of entanglement appear in even the 
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most elementary measurements (e.g. interference of Ry- 
dberg series in photoabsorption spectra). A quantitative 
measurement of entanglement, that would more closely 
reflect the evolution of the linear entropy, can be set up 
by combining the methods empl oyed for the detection of 
interfering Rydberg wavepackets |22| . These methods are 
based on the use of several laser pulses with well-defined 
phase relations to monitor the interferences appearing in 
the population of predefined Rydberg states. 

V. CONCLUSIONS 

We have investigated the entanglement production in 
Rydberg molecules and contrasted the results with the 
underlying classical dynamics, which is known to play 
a role in the understanding of observable spectroscopic 
effects and in the interpretation of the energy lev- 
els statistics We have first determined the aver- 
age linear entropy of eigenstates corresponding to dif- 



ferent collision strengths (quantum phase shifts or clas- 
sical kicks) and then studied the generation of entangle- 
ment from initially unentangled states. We have seen 
that the quantum/classical correspondence on the level 
of entanglement production is relevant not on the scale 
of the global classical dynamical regime, but rather on 
the specific classical features that quantum-mechanically 
translate into superpositions. In Rydberg molecules, it 
is a high rate of inelastic scattering for an initial clas- 
sical distribution that corresponds in the quantum do- 
main to states displaying a high degree of entanglement. 
The relation between global chaotic or regular behav- 
ior on the one hand and these specific classical features 
that will be translated quantum mechanically as entan- 
glement production on the other hand appears to highly 
depend on the individual system under investigation. We 
therefore conclude it seems unlikely that the generation 
of entanglement could be employed as a reliable signature 
of chaos for an arbitrary dynamical system. 
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